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ABSTRACT 

We have studied Very Long Baseline Array (VLBA) polarimetric observations of 8 
sources including quasars and BL Lacs at 12, 15, 22, 24 and 43 GHz and high frequency 
rotation measure {RM) maps are presented. We find typical values for the RM in the 
VLBI core of several thousand rad/m^, which are higher than values in the literature 
at lower frequencies. Assuming a dependence of the form RM oc z^^, we obtain an 
average value of a = 3.6 ± 1.3, which is larger than that expected by theoretical 
considerations. Rotation measures are detected in the jet of only two sources and 
we find that only 0906+430 (and possibly 1633+382) show indications of a robust 
gradient. We discuss the Faraday-corrected polarization properties of the sources. Our 
interpretation supports the presence of helical magnetic fields with new, unresolved, 
components affecting the intrinsic direction of polarization close to the base of the jet 
of some objects. 

Key words: Galaxies: active. 



1 INTRODUCTION 

The emission associated with jets in active galactic nu- 
clei (AGN) is due to synchrotron radiation over a wide 
range of the electromagnetic spectrum, including radio fre- 
quencies. The emission is generally highly linearly polar- 
ized . more than 50% in so me optically thin jet regions (see 
e.g. i Lister HomanlliooHh . clearly indicating the presence 
of highly ordered magnetic fields. In several cases it has 
been found that this polarization is perpe ndicular to the 
direction of the jet taabuz da et al.l l2Q0d and references 
there i n) or has a spine/sh eath structure (|Attridge et al.l 
I1999I : IPushkarev et al.l |2Q05|): see also numerous examples 
in the MOJAVE database (jLister Homanll2005l ). 

Although this can be interpreted as a series of shocks 
in the jet enhancing the perpendicular component of the 
magnetic field, and /or the inter action of the jet with the 
surrounding media (lLaindll980h . the presence of polariza- 
tion either parallel or perpendicular to the jet can also be 
easily explained with a helical magnetic field. There is evi- 
dence that helical magn etic fields seem to be present in, at 
least, a fraction of AGN. iGabuzdal (|2006l ) summarizes differ- 
ent observational polarization tests in order to distinguish 
between helical magnetic fields and alternative explanations 
concluding that there are a number of sources for which rea- 
sonable evidence for helical jet magnetic fields is available. 

One way to study the structure of the intrinsic mag- 
netic field is the analysis of the rotation measure (RM) dis- 
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tribution across the source. The polarization angle x rotates 
following the relation x Xo + RMX^, where xo is the in- 
trinsic polarization angle, A is the wavelength and RM is 
the rotation measure, given by 



RM oc 



B-dL 



(1) 



where Ue is the electron density and B-dl is the magnetic 
field along the line of sight. 

The first implication is that the intrinsic polarization 
(and hence the inferred direction of the magnetic field giv- 
ing rise to it) will be altered by the effects of the Faraday 
rotation. As the RM depends on the component parallel to 
the line of sight of the magnetic field, we will have, in gen- 
eral, a variable alteration of the polarization direction across 
the source. Thus, if we want to properly understand the in- 
trinsic properties of the magnetic field, we need to map the 
RM across the source in order to adequately subtract its 
effects in the different regions. 

If we have a helical magnetic field, its toroidal compo- 
nent will produce a rotati on measure chan ge across the jet, 
giving rise to a gradient (|Blandfordl[l993l ). In the simplest 
case, with the jet perpendicular to the line of sight in the 
observer's rest frame, we would observe an antisymmetric 
distribution of the rotation measure, with positive values on 
one side of the jet, negative values on the other and a null 
RM in the centre. As the viewing angle decreases, there will 
be an offset of the absolute value of the RM values and, 
if the viewing angle turns out to be smaller than the pitch 
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angle, we will observe only positive values for the rotation 

measure (lUchida et al.ll2QQ4l : IXsada et al.ll2002l). 

Such RMs have already been detected: *Zavala Tavloi 

20031 . [2004) studied RMs in a sample of AGNs. Asa da et al 



2002 



20081 ) found a time-variable RM gradient over 



more than lOOpc along the jet in 3C273, the origin being 
quite probably t he sh eath around the ultra-relativist ic jet. 
iGabuzda et "aU (|2004) also found rotation measure gradi- 
ents ranging from negative to positive values across the jet 
in 0745-h241. However, the stability o f the RM structures 
is still unclear: IZavala Tavlorl (|200lh found time variabil - 
ity of the RM in 3C273 and 3C279 but lGomez et al.1 (l201lh 
found that, although there were changes in the linear polar- 
ization of 3C120, the underlying rotation measure remained 
unaltered along 2 to 5 mas from the core over several years. 

Previous results are based on observations obtained 
with the Very Long Baseline Array (VLBA) using, with the 
exception of few cases, low frequency bands (typically in 
the range of 8-15 GHz) based on a compromi s e bet ween 
resolution and sensitivity, and only iGomez et al.l (|201lh and 
lAttridge et al.l (|2005l ) have obtained results for frequencies 
as high as 86 GHz for the radio galaxy 3C120 and 3C273, 
respectively. Thus, we seek to obtain more measurements at 
higher frequencies so that we can study the behavior and 
structure of the rotation measure closer to the base of the 
jet in a larger number of cases. We can then compare our 
results with the ones obtained at lower frequencies. 

In this paper we study the rotation measure of 8 sources 
observed with the VLBA from 12 to 43 GHz. In section 2 we 
summarize the observations and data reduction. In section 
3 we present our results. In section 4 we discuss our findings 
and we present our conclusions in section 5. An analysis of 
the robustness of rotation measure gradients is presented in 
the appendix. 



2 OBSERVATIONS AND DATA REDUCTION 

Polarization observations of 6 radio-quasars (0133+476, 
0420-014, 0745+241, 0906+430, 1633+382 and 1954+513) 
and 2 BL Lac objects (0256+075 and 1823+568) were car- 
ried out in a 24 hour session on November 2, 2008 using the 
VLBA. The frequencies selected were 12.039, 15.383, 21.775, 
23.998 and 43.135 GHz, each with 2 IPs, with a bandwidth 
of 8MHz and a bitrate of 128MB/s. The sources were ob- 
served in a "snapshot" mode, with 8-10 several minute scans 
for each frequency and object spread in time, so that the re- 
sulting UV coverage was quite uniform. The data reduction 
and imaging was done with the NRAO Astronomical Im- 
age Processing System (AIPS) using standard techniques. 
The reference antenna used was Los Alamos. Simultaneous 
solutions for the instrumental polarizations and source po- 
larizations for the compact source 1954+513 were derived 
using the AIPS task LPCAL. 

We calibrated the electric vector position angles (EV- 
PAs) using the VLBA D-Terms (Gomez et. al, 1992). Com- 
parison of the D-terms against a set of tabulated values 
(previously calibrated by other means) proves to be a reli- 
able method for calibrating the absolute L-R phase offset 
in VLBA observations. We were able to obtain two sets of 
calibrated data with either the same or nearly the same fre- 
quencies (S. P. O 'Sullivan, private communication; A. Re- 



ichstein, private communication) so that it was possible for 
us to apply this method to our data. 

To ensure the reliability of this method, we performed 
a series of extra checks. First, we compared the indepen- 
dent results given by the two comparison sets of D-terms 
for consistency. Second, we compared our new images with 
our previous images (Algaba 2010). Third, for some frequen- 
cies, we were also able to compare the resulting polarization 
angles with images from the MO JAVeQ, University of Michi- 
gan Radio Astronomy Observatory Databas43 and the 7-mm 
monitoring by Boston University Blazar Group databases. 
Fourth, we obtained the EVPAs of several sources in the 
optically thin jet regions and checked for consistency with 
modest Faraday rotation. We estimate the overall error in 
the EVPA calibration to be about 4° for all frequencies. 

After the initial construction of total intensity maps, 
we made new versions at 12, 15, 22, 24 and 43 GHz us- 
ing the best calibrated data but convolving all frequencies 
with a circular version of the beam obtained for the 15 GHz 
map. This beam was chosen as a compromise between the 
lower (12 GHz) and higher (43 GHz) resolution maps, super- 
resolving 12 GHz only by a modest amount. We then made 
the maps for the Stokes parameters Q and U that we used to 
construct the polarized flux (p = -\- W^) and the polar- 

ization angle (x = ^arctan^), with their respective noise 
maps, using the AIPS task COMB. As information about 
the absolute position is lost during the calibratio n process, 
and given the core shift at different frequencies ([Lobanovl 
1 19981 ). we had first to properly align the different maps at 
various bands. For this, w e used the program developed by 
'Cro ke Gabuzdal (|2008l ) based on the cross-correlation of 
optically thin regions of the jet. 

For the construction of the rotation measure maps, as 
the current version of the AIPS task RM is limited by a 
maximum of 4 input frequencies, we used a modified version 
obtained from R. Zavala that can construct rotation measure 
maps using up to 10 frequencies. The distribution of RM was 
obtained in the regions where the error in the input maps 
was smaller than 15°. We found that the results are very 
similar compared to when we blank those regions with the 
error on the RM exceeding 100 rad/m^, except the blanking 
edges appears smoother in the former case. In order to check 
the goodness of the A^-fit we used the AIPS task RMCUB, 
also included in R. Zavala's package, which allows us to plot 
a grid of RM fits and obtain fit parameters, such as RM fit 
error, x^, correlation coefficient and the Q-value. 

For 0420-014, 0745+241, 1823+568 and 1954+513 we 
used all 5 av ailable frequencies. T he case of 0133+476 is very 
complex fsee lAlgaba et al]|2012l ). including possible integra- 
tion of different polarization components or frequency vari- 
ability of the rotation measure, and thus we decided to only 
fit our three highest frequencies. In the cases of 0256+075 
and 0906+430, the polarization structure becomes more 
complicated when reaching the highest frequency in our ob- 
servation. Consequently, the 43 GHz EVPA value signifi- 
cantly departs from the trend seen in other frequencies. This 
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indicates that we are either probing a different region or that 
we are integrating over different non-resolved components. 
Hence, we consider these values not to be representative of 
the region we study in the rest of the frequencies and so 
we decided to not include 43 GHz in our RM maps. For 
1633+382 the 43 GHz angles have larger errors, thus we 
exclude this frequency from our study. Frequencies used in 
each source for the calculation of RM maps are summarized 
in Table 1. 

There are a number of factors which might alter the ori- 
entation of the polarization angles in our maps. For example, 
a transition in the optical depth of the source can induce a 
90° change in the polarization angle. Also, 180° ambiguities 
ar e inherent to the pola rization angles. We investigated these 
m lAlgaba et all (|2012l ). where we analyzed the total and 
polarized intensity properties of the sources, finding that all 
sources were either only optically thin or only optically thick 
in the regions where we find RM. Hence, we are confident 
that no 90°jump due to a thin-thick transition is invalidat- 
ing our results. In the same way, we studied the possibility 
of 180° ambiguities in the polarization angles in the region 
close to the core and these have been taken into account in 
the calculations. 

In order to double check our results, various RM maps 
were created for each source using different blanking meth- 
ods and/or a different set of frequencies, comparing the fits 
both by eye and checking the goodness of the fit using the 
task RMCUB. We noticed that, for sources like 0256+075 
this implied the addition of noisy areas in the RM map and, 
in some cases as in 0133+476, the emergence of patchy areas. 

In cases such as 0256+075 and 0906+430, including 43 
GHz resulted the fit not to pass through all the data and 
tilted the slope to a wrong value. As we excluded 43 GHz 
in these, the goodness of the fits improved from R^ ^0.35 
and 0.86 to R^ -0.95 and 0.96 for 0256+075 and 0906+430 
respectively. For other sources, such as 1633+382, the RM 
map did not change significantly, but due to the error asso- 
ciated with the additional frequency, the fit turned out to 
be worse. 

We also checked the consistency of the RM maps when 
a shift on any of the input maps was applied. For this, 
we considered input maps with and without a core shift 
(iLobanovll 19981 ). as well as horizontal or vertical shifts up to 
200/ias, and we found no significant difference in the output 
RM maps. 



3 RESULTS 

Rotation measure maps derived for all 8 sources studied here 
are shown in the top panels of Figure 1. Here, the contours 
correspond to the total intensity of the convolved version of 
the map for 15 GHz, and the color levels indicate the ro- 
tation measure found. Peak and bottom intensity contours, 
corresponding to 3xRMS, are given in Table 1, together 
with the frequencies used for the Faraday rotation A^— fit, 
and the rotation measure scale for all sources. In all maps, 
contours increase in steps of a/6, the color code for the ro- 
tation measure is shown in the top, and the beam in the 
bottom left corner of the image. 

Rotation measure is typically found in the VLBA core 
region; that is, close to the base of the jet. In 0745+241 



some RM extending a bit further from the core is found, 
and only in 0906+430 the RM extends continuously from 
the core up to a distance of 2 mas into the jet. In 1633+382 
and 1954+513 we also find RM in the knots located at 4 and 
3.5 mas from the core, respectively. Typical RM values are 
of the order of some krad/m^ with 1633+382 and 1954+513 
showing the most extreme values. 

In order to investigate the observed RMs we made slices 
to analyze the rotation measure across the jet. Plots of the 
rotation measure across the jet are shown below the RM 
maps in Figure 1, indicating the region taken for the slice 
with a thick black line in the map (note that for 0420-014 
and 0745+241 no slices are shown, as these sources do not 
present any indications of variations of the RM across the 
jet). In all cases, the direction for the slice is taken to be 
from left to right and, in case of confusion (1633+382, for 
example), from top to bottom. The horizontal axis of the 
rotation measure slice plot has been zeroed at the beginning 
of the slice. A line in the bottom of the slice indicates the size 
of the beam. Note that, for the case of 0133+476, the slice 
is smaller than the beam size. We discuss the significance of 
these gradients following the criteria given in the Appendix. 

0133+476 is a quasar showing the RM in the bound- 
aries of the VLBI core. The values obtained for the RM span 
from -4000 to 1000 rad/m^. Although the slice taken shows 
a very clear gradient with a large SNR (cr — 15), the region 
where a possible gradient is seen is very limited, less than 
one beam width, making this very ambiguous. According to 
the simulations performed in the Appendix, and together 
with the fact that it is located in the core, this is not a 
reliable gradient. 

0256+075 shows a gradient in the core region spanning 
a bit more than a single beam size and forming an angle with 
the direction of the jet. The values for the RM found in this 
source vary monotonically from 2000 rad/m^ in the upper 
region of the core to about 1000 rad/m^ as we approach 
the upstream region of the jet. However, this gradient is too 
small both in size and slope to be considered reliable. 

0420—014 shows no indications of variations transverse 
to the jet but does show the RM monotonically decreas- 
ing from positive +1000rad/m^) to negative {^^ —1000 
rad/m^) values as we move further from the central region. 
0745+241 shows a behavior similar to these of 0420-014 
with higher RM values (^^3000 rad/m^) closer to the core 
and a monotonic decrease to ^1000 rad/m^ at a distance of 
^1 mas. 

0906+430 has an almost null RM in the core, whereas 
it takes positive values (+2000 rad/m^) on one side of the 
jet, decreasing gradually and changing sign down to -2000 
rad/m^ on the other side of the jet. This is the source which 
shows the clearest RM gradient spanning two beam widths 
along the jet. The slope is about 3 times the typical er- 
ror and its profile is very smooth, in agreement with all 
the requirements discussed. Also, the gradient is clearly vis- 
ible at around one beam size away from the core, in agree- 
ment with simulations performed bv lBroderick McKinnevI 
(l201Ql ). We consider this RM gradient to be robust. 

1633+382 Shows a very high (RM=22 krad/m^), fea- 
tureless RM in the core, but displays an interesting gradient 
in the knot located at about 4 mas from the core. We discuss 
reliability of this gradient in section 4.3.2. We note that in 
Figure Id we have set the RM range to enhance the RM 
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(a) 0133+476 (b) 0256+075 (c) 0906+430 



-0.5 0-0 0^5 1^ 0.6 0^8 1^0 1^2 t^^ -15.0 -14.5 -14.0 -13.5 -13.0 
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(d) 1633+382 (e) 1823+568 (f) 1954+513 



Figure 1. Rotation measure maps. For each of the sources, top: RM map; bottom: shce of the RM along the thick black line taken on 
the RM map. Color figures can be found in the online version. 
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Table 1. Parameters used for the maps. 



Source I Peak Bottom Contour Polarization Sticks Frequencies used for RM RM range 

(Jy/beam) (mJy/beam) (mJy/beam)/mas (GHz) (rad/m^) 

[1] [2] [3] [4] [5] [6] 



0133+476 
0256+075 
0420-014 
0745+241 
0906+430 
1633+382 
1823+568 
1954+513 



3.19 
0.32 
3.72 
1.39 
0.69 
2.18 
1.18 
1.21 



3.2 
1.0 
5.6 
1.4 
1.4 
2.2 
2.4 
1.2 



50 
50 
200 
40 
10 
20 
100 
40 



22, 24, 43 
12, 15, 22, 24 
12, 15, 22, 24, 43 
12, 15, 22, 24, 43 
12, 15, 22, 24 
12, 15, 22, 24 
12, 15, 22, 24, 43 
12, 15, 22, 24, 43 



-4000 - 1000 
1000 - 2000 

-1000 - 1000 
500 - 3000 

-2000 - 2000 
-500 - 1000 
500 - 1500 
-15000 - -13000 



features in the knot for clarity in the discussion. As high RM 
values are not included in the scale, the map is saturated for 
the values occurring in the core. 

1823+568 shows a very flat rotation measure in the 
core, with values of ^1200 rad/m^. As we move away from 
the VLBI core some indications of a transverse gradient tak- 
ing only positive values arise. This is in general agreement 
with the idea that finite-beam effects blend different regions 
and make the detection of reliable features of the RM in the 
core impossible. 

1954+513 shows indications of a gradient across the 
core region that seem to be faded out (or inverted) as we 
move less than 1 beam si ze away. This case is similar t o 
one of the simulations of iBroderick McKinnevI (|2010l V 
where spurious rotation measures and sign changes occur, 
and thus, we do not consider this gradient reliable. 

In general, we find rotation measures to be on the order 
of several krad/m^, the core of 1633+382 being the most ex- 
treme case. Variations of the RM in the optically thick core 
are, in general, not reliable and give rise t o unrealistic gra- 
dients (IBroderick k McKinnevI l2010l ) and iTavlor k Zavalal 
Thus, such RM variations will not be discussed. 



4 DISCUSSION 

4.1 Comparison with previous results 

The source 2 56+075 has been previously studied in 
iMahmud et al.l (|2009al ), where they claimed to have found 
an RM gradient in the VLB A core, with a gradient reversal 
th at they explained u s ing m agne tic tower models pr oposed 
bv lLvnden-Bell et al.1 (|l994h and iLvnden-Belll (Il996l ). How- 
ever, their gradient was found very close to the core and with 
a width not significantly exceeding the beam size. Thus, it 
is likely that the gradients they observe are also due to RM 
fluctuations in the core. 

0420-014 was one of the sources for which an RM anal- 
ysis was done usin the MO J AVE p rogram. Almost no RM 
was found by Ho vatta et al.l (|2012l ) for this source in Oc- 
tober 2006, with only small patchy traces within one beam 
size from the core. The difference in the RM detection could 
be due to the increase of both polarization and the degree 
of polarization by a factor of 2 in the epoch we study here, 
allo wing us to obtain a b etter SNR with a similar sensitivity. 

IVenturi k Tavlorl (|l999h performed observations of 
0906+430 in November 1996 and also found clear RM over a 
long distance from the core in this source. The rotation mea- 
sure they found in the core is on the order of 100 rad/m^, 
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about one order of magnitude smaller than the results pre- 
sented here. No R M of this source is sh own in the Faraday 
rotation studies of iHovatta et alJ (|2Q12l ). 

Rotation measure in the sources 133+476 and 
1823+ 568, and 0745+241 was studied by iHovatta et "aD 
(|2012l ) and lGabuzda et alJ (|2004l ) respectively. The RM fea- 
tures and sign are different when compared with our work. 
Although this could be due to some time variation, given 
their location (optically thick core) and extension (generally, 
less than 1.5 beam sizes), it is also likely that they might 
have a non-negligible contribution from spurious RM due 
to noise in the data. 

4.2 Increase of Core RM with Frequency 

As we observe at higher frequencies, the core RM is expected 
to increase. This is due to the apparent change in location of 
the central engine with frequency (i.e, core shift) because of 
optical depth effects. Thus, observations of the VLBI core at 
higher frequencies probe regions closer to the central engine, 
where both magnetic field and particle densities increase. 
Hence, the core RM is expected to increase with fre quency 
following the formula RM oc ly"" (jjorstad et al.ll2007f ). The- 
oretical estimations provide a = 2 under the assumption of 
a toroidally dominated magnetic field and electron density 
in equipartition scaling as B oc and Ue oc , respec- 
tively, where d is the distance from the central engine, and 
the outflow is a spherical or conical wind. This value is also 
in general agreement with estimatio ns for a on oth er similar 
sources bv lJorstad et^ (2007) and[o: Sullivan Gabuzda 

liooi). 

We have compared our measurements with previous 
core RM results in the literature at other fre quencies. For 
our d ata, we used the core RM obtained in lAlgaba et al.l 
(|2012l ). For data from the literature, we used the author's 
fitted data for the core. If unavailable, we estimated the core 
RM base d on their RM map. We note that, for the case of 
0420-014, IHovatta et"all (l2012l ) were not able to find RM in 
the VLBI core and thus we have used the value they found 
upstream in the jet as the lower limit. To our knowledge, 
there is no previous VLBI RM data for 1954+513. 

Our results are summarized in Table 2, where column 
1 indicates the source; columns 2, 3 and 4 the reference of 
the lower frequency RM data, the frequency used in the 
literature and the core RM they found for that frequency, 
respectively; columns 5 and 6 the lowest frequency we use 
from our data to determine the RM and its value; and col- 
umn 7, the value for a we obtain when we compare previous 
data with ours at higher frequency. 

We find values obtained are around a = 3, with a lower 
limit being a < 1.2 for 0420-014 and the highest value a = 
11.2 for 1633+382. If we exclude these two values, we find 
the average (a) = 3.6+1.3. Except for the case of 1633+382, 
with exceptio nally high a, these values are similar to the 
ones found bv 10 'Sullivan Gabuzdal (l2009l ). ranging from 
0.9 < g < 3.8, but a bit larger than the ones in ljorstad et al.l 
(|2007l ). The average value found here is not compatible with 
the theoretical estimation a = 2 and, indeed, only the value 
for 0133+476 matches this value to within la. Two possible 
explanations for this include time variations or convolution 
effects. 

Our observations were taken roughly two years later 



than several of the lower-frequency observations we are com- 
paring with. During this time, several changes may occur in 
the innermost regi ons of the je t causing a change of the core 
RM. Indeed, Zavala Tavlorl dioOl. ) found changes in the 
core RM in 3C273 and 3C279 during a period of about 1.5 
years. In our case, this is particularly true for 0745+241, 
0906+430 or 1823+568, sources that show a change in the 
sign of the RM for the two epochs. Such RM sign changes 
have been detected before for 0133+382, 0745+241 and 
1633+382 when com pared with another set of observations 
(jAkaba et al.ll2012l ). 

A different reason for our relatively high a value could 
be that convolved RMs are, in general, smaller than the true 
values and a function of the beam size (iMurphv Gabuzdal 
I2OI2I ). As we observe at higher frequencies, the beam size 
becomes smaller and, thus, this effect becomes less impor- 
tant. In practice, this implies that the RM observed at lower 
frequencies was underestimated (or, at least, more so than 
the RM values studied here), thus producing an artificial 
increase of the a parameter. However, the range of observed 
frequencies is relatively small and thus the effects due to de- 
convolution on the RM are expected to be similarly small. 

Another possibility is that the assumptions for the theo- 
retical derivatio n of a = 2 are not ade quate for these sources. 
As discussed in lAlgaba et al ] (120121 ). it is a reasonable as- 
sumption to consider that these are in equipartition. How- 
ever, it is not yet clear how magnetic and electron densities 
actually scale with distance from the central engine. If these 
scale in a different way, this could lead to a theoretical pre- 
diction for a that is closer to our measured value. For ex- 
ample, if the electron density decays as Ue oc d~^, with the 
rest of the parameters left unchanged, then the theoretical 
estimate would be a = 3. Further investigation of RMs at 
different wavelengths should be performed in order to clarify 
this. 

4.3 Rotation Measure Gradients 

Although rotation measures have been studied for several 
decades, it is still not clear when observed variations and 
gradients of RM are reliable and only re cent studies and 
simulations have dealt with this question. iTavlor Zavalal 
(Hoio) introduced a series of criteria to dis c ern w hen a gra- 
dient was reliable, Broderick McKinnevI (|2010h discussed 
jet RMs from theoretical ma gneto-hydrodynamic simula- 
tions and IHovatta et al.l (|2012l ) simulated errors of Faraday 
rotation for their observations. We discuss these and perform 
similar simulations in the appendix. 

In general, the RMs we find are located in the 
core of the sources. In this region, the RM structure 
i s less reliable and is subjec t to several blending effects 
(jBroderick McKinnevI |20 id ). Thus, we will not consider 
variations of core RM. We have found indications of a reli- 
able RM across the jet in two of the sources studied here: 
0906+430 and possibly 1633+382. 

The quasar 0906+430 displays a constant RM gradient 
that spans along the jet up to a distance of 2 mas. At a 
redshift of z — 0.67 and its cosmology corrected scale 6.831 
pc/mas, this implies a rotation measure gradient extend- 
ing at least '^14pc. The amount of rotation of the polariza- 
tion angles that we find is greater than 90° (|Algaba et al.l 
l2012h and the fractional polarization, although relatively 
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Table 2. RM increase with Frequency 



Source 


Reference 


1^1 


Core RM{ui) 


1^2 


Core RM{u2) 


a 






GHz 


rad/m^ 


GHz 


rad/m^ 




[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


0133+476 


H12 


8 


-216 


22 


-2500 


2.4 ±0.4 


0256+075 


M09 


4 


50 


15 


1530 


2.6 ±0.5 


0420-014 


H12 


8 


> 548 


12 


880 


< 1.2 


0745+241 


G04 


5 


-120 


12 


2550 


3.5 ±0.5 


0906+430 


V99 


4.8 


100 


12 


-3200 


3.8 ±0.5 


1633+382 


H12 


8 


-235 


12 


22040 


11.2 ±0.4 


1823+568 


H12 


8 


-121 


12 


1250 


5.8±0.4 


1954+513 








12 


-13800 





G04: Gabuzd a et al.l (l2004l ): V99 : IVenturi Taylor i (il999l ): 
H12: Hovatta et al.l (|2012l ): M09: iMahmud et al.l (|2009af ) 



weak (< 1% at 12 GHz), is smoothly increasing with fre- 
quency. Hence, we can estimate that the contribution of in- 
ternal Faraday rotation, if any, is negligible. In turn, this 
indicates that the screen producing the Faraday rotation 
gradient must be external. 

The quasar 1633+382 shows a rotation measure at 
about 4 mas from the core. If we take into account the er- 
rors, the RM runs from around a thousand rad/m^ to zero 
(or possibly negative values). Unlike the case of 0906+430, 
this RM is located in a knot of the source, where both total 
intensity and degree of polarization are enhanced. The range 
of the RM indicates that, as in 0906+430, it is dominated 
by external Faraday rotation. 

The case of 1633+382 is indeed close to the limits of 
what we could consider a reliable gradient: it spans over 
about 1.5 beams across the jet with a significance a 2. As 
we state in the appendix, there is a chance of about ^ 1% 
for this gradient to be a false detection. It is possible to 
have spurious gradients of up to ^230 rad/m^/beam over 
1.5 beams (c.f. Figure A4b). This would lead to up to a 
contribution of ^345 rad/m^ due to the noise, which is about 
25% of the gradient observed in this source. Therefore, the 
existence of a physical gradient is possible, although we shall 
discuss it with caution. 

RM is al so found in th e jet o f this source in the Fara- 
day studies of lHovatta et al.l (|2012h . Typical values they find 
in the region we are studying here are RM ^ 200 rad/m^, 
which is ab out five tim es smaller than the values we find. In 
their studv, [Hovatta e t al. (20li) do not find strong indica- 
tions of a gradient. If the RM gradient we detect is real, we 
suggest that it might be due to evolution of the region (see 
below) . 

4.4 Polarization Structure of the Jets 

Once we have obtained RM maps, it is possible to correct 
for Faraday effects on these regions. We show Faraday cor- 
rected electric field maps in Figures 2-4. Here, contours in- 
dicate total intensity at 12 GHz as in Figure 1 and direc- 
tion and longitude of the sticks indicate the direction of the 
Faraday corrected E-field direction and polarization inten- 
sity respectively. Note that we do not show the E-field for 
these regions where we did not find RM ^ as we are unable to 
correct it for Faraday rotation. A longitude of 1 mas in the 
polarization sticks corresponds to the number of mJy/beam 
indicated in column 4 of Table 1. 



At a first glance, it seems that variation in the direction 
of the electric field in the core region is smaller and smoother 
when corrected for Faraday effects. This is particularly true 
for the sources 0256+075, 0420-014 and 1823+568. How- 
ever, in sources such as 0133+476, 0745+241 or 0906+430 
there appear to be two differentiated components in the di- 
rection of the E-field: one in the VLBI core and one up- 
stream in the jet. We suggest that this second component is 
caused by an unresolved feature in the jet. 

In several sources, the electric field in the core seems to 
be either perpendicular (0745+241, 1954+513) or parallel 
(0133+476, 0256+075, 0420-014, 0906+430) to the direc- 
tion of the jet. This bi modal configuration h as been previ- 
ously observed (see e.g. iGabuzda et ahlliood and references 
therein) and can be interpreted as reflecting the presence 
of helical magnetic fields associated with the jets. We ob- 
serve polarized emission from both the front and the rear 
of the emitting regions but the polarization angles will have 
a 180° offset and, when both contributions are vectorially 
added, the result will depend on the pitch angle of the jet. 
With a small pitch angle, the azimuthal component of the 
magnetic field dominates and the vectorial addition will pro- 
duce polarization perpendicular to the jet; whereas with a 
large pitch angle, the toroidal compone nt dominates and t he 
polarization will be parallel to the jet (|Asada et al.ll2003 v 

In only two sources (1633+382 and 1823+568) there 
is an evident misalignment between the E-field and jet 
directions. One of the explanations for this, as above, is 
the presence of an unresolved (in this case even in polar- 
ization) component. Other possibilities include opacity ef- 
fects on the magnetic field direction or differential Doppler 
boosting, which would twist the electric vector pattern 
([Roberts Wardlell2012l ). Oblique shocks with arbitrary an- 
gles can also cause arbitrary observed polarization angles. 

Due to sensitivity limits, we have been able to detect 
significant polarization at high frequencies and RM in the 
jets of only two sources: 0906+430 and 1633+382. We study 
these two cases below. 



44.1 0906+430 

Figure 3 shows the polarization structure of this source. The 
top panel shows the polarization map as detailed above. The 
bottom panel shows the degree of polarization m of the map 
at 12 GHz across the width of the jet in the direction of 
the RM slice shown in Figure 1. The axis has been zeroed 



8 



J. C. Algaba 




a 




-O 








(a) 0133+476 



(b) 0256+075 



(c) 0420-014 



n \ r 






MilliARC SEC 



MilliARC SEC 



MilliARC SEC 



(d) 0745+241 



(e) 1823+568 



(f) 1954+513 



Figure 2. Faraday corrected magnetic field maps for sources with RM found only in the core. Contours indicate total intensity (as in 
Fig. 1), whereas sticks length and direction indicate strength and direction of E-field (see Table 1, column 4). 



at the centre of the jet. In the degree of polarization slice, 
areas shaded in grey indicate the regions with polarized flux 
density below 3xRMS in the polarization images. 

There seem to be two regions with different electric field 
directions along the jet. Close to the VLBI core, the E-field 
appears to be forming an angle with the direction of the 
jet, whereas further downstream it becomes perpendicular 
to the direction of the jet, as in the core. As discussed above, 
this might be caused by an unresolved component, opacity 
or relativist ic effects. If we investigate the 43 GHz intensity 
map, we can identify a resolved component located at 1 
mas from the core, and hence we think this might make 
a significant contribution to the change of the polarization 
direction. 

There are previous studies (see e.g. iMarscher et. al.l 

l201Ql ) where a rotation of the optical polarization angle has 
been seen, presumably due to a structure following a spiral 
path through a helical magnetic field where the flow accel- 
erates. One of the possibilities is that we are observing a 
"snapshot" of the radio counterpart here, with this struc- 
ture rotating the polarization angle or the region around 
it. 

When we study a cut of the degree of polarization m 
along the direction of the gradient (i.e, across the jet), we 
find a minimum in the centre of the jet and indications of an 
increase as we move toward the edges, although the errors 
become rapidly large and, after 0.5 mas from the jet cen- 



tre, the errors dominate and the polarization drops below 
3xRMS. In the region of the RM gradient, even if we take 
into account the errors, m seems to have a concave shape 
across the jet. 

Overall, one possible way to explain the structure of 
0906+430 could be as follows: a new component (only re- 
solved at 43 GHz in total intensity) is emerging from the 
VLBI core and following the path along the intrinsic helical 
magnetic field, thus altering the direction of the observed 
polarization closer than 1 mas to the core. Further away, we 
find three aspects (RM gradient, spine structure of the po- 
larization and shape of m across the jet) that point toward 
ajhelical shape of the magnetic field. We note (jContopoulosI 
l2009l ) that a spine-sheath polarization structure is not a nec- 
essary condition to pinpoint a helical magnetic field, as the 
observed structure depends on parameters such as viewing 
and/or pitch angles. 



4.4.2 1633+382 

Figure 4 shows the polarization structure of this source. The 
top panel shows the polarization map; the middle panel, the 
degree of polarization m of the map across the width of the 
jet in the direction of the RM slice zeroed at the centre 
of the jet. The bottom panel shows m along the jet zeroed 
at the core. As before, in the degree of polarization slices. 
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Figure 3. Polarization of 0906+430. Top: polarization map. Bot- 
tom: degree of polarization across the jet following the direction of 
the RM slice. Grey indicates the areas with polarized flux density 
below 3xRMS in the polarization images. 



areas shaded in grey indicate the regions with polarized flux 
density below 3xRMS in the polarization images. 

In the region of the jet where we have been able to 
correct for Faraday rotation, the E-field seems to form an 
angle with the direction of the jet. When we take a slice 
perpendicular to the direction of the jet to study the polar- 
ization structure of this feature, we find hints of a complex 
W-shape, although, given the errors both in total and frac- 
tional polarization, the only significant feature is a peak of 
the fractional polarization in this region. 

A smooth sheath around the jet cannot produce a region 
of enhanced polarization. Hence, we speculate, based on the 
total intensity and polarization maps, that the reason for 
this shape is a component arising from the polarization of 
the knot itself. Indeed, the local peak of the degree of polar- 
ization is coincident with the intensity peak of the knot, and 
both shapes also agree within the errors. We note that this 
is not unique, and has also been observed in other sources 



60 



40 - 



20 




60 



40 - 



20 - 



-1 



/ 



12 3 4 

Section Along the Jet [mas] 



Figure 4. Polarization of 1633+382. Top: polarization map. Mid- 
dle: degree of polarization across the jet following the direction of 
the RM slice. Bottom: degree of polarization along the jet. Grey 
indicates these areas with polarized flux density below 3xRMS 
in the polarization images. 



10 J. C. Algaba 



l|Algaba et all 120131 '). We also note that magnetic fields do 
not need to be perpendicular to the direction of the jet to 
be caused by a shock: oblique shocks with arbitrary angles, 
leading to different directions of the magnetic field are also 
possible. 

In order to check if such a knot is compatible with the 
observations, we have obtained a slice of m along the jet. We 
observe a very low m near the core, but a peak of m ~ 20% 
at a distance ^ 3.5 mas (corresponding to ^700 pc, depro- 
jected), coincident with the location of the knot. According 
to iLiul (|201Q| ) and MOJAVE data, this component has been 
moving at a velocity of ^ 0.32 mas yr~^, gradually deceler- 
ating until reaching its current position around 2003, where 
it seems to remain steady. 

One possibility is that this is a component that has been 
gradually decelerated by interaction with the surrounding 
media and has developed a shock front with compressed 
magnetic field. Given the degree of polarization, we can 
conclude that the underlying magnetic field is highly or- 
dered. Compression across the jet will produce an enhance- 
ment of the t ransv erse component of the magnetic field. 
lAkaba et all (|2012l ) found that, for this source, the mag- 
netic field at Ipc from the central engine is B{lpc) = 0.7G. 
If we assume B oc which is in agreement with equipar- 
tition arguments and can be used to describe the observed 
synchrotron emission in compact regions of VLBI jets, we 
find that the corresponding magnetic field strength in the 
region of the knot is around B{lpc) ^ ImG. This is a rela- 
tively high magnetic field strength but is in agreement with 
the one found in other sources at similar scales (see e.g. 
lOwen et al.lll989l ). 

There is another possibility that could also give rise to 
an intrinsic Faraday-corrected polarization forming an ap- 
parent angle with the jet in this region. If we analyze the 
structure of this source at higher frequencies, a component 
close to the core ('-^ 0.5 mas) is observed in the north-west di- 
rection (see e.g. the 7-mm monitoring by Boston University 
Blazar Group). This could be an indication that the direc- 
tion of the jet has actually changed and new components are 
now being ejected in a direction that is tilted with respect 
to the previous one. Alternatively, the co mponent l ocated 
at 0.5 mas may be a standing component (iLiulbOlol ). thus 
indicating the truly intrinsic direction of the jet and that the 
actual jet is curved. If this is the case, polarizations both in 
the core and in the knot might be showing the true direction 
of the jet in these regions (i.e., the E-field would be parallel 
to the direction of the jet). 

We have obtained new slices of m assuming that the 
local direction of the jet is the one shown by the E-field 
vectors and we find that the previous discussion does not 
significantly change: namely, the W-shape across, and the 
peak along, the local direction of the jet in the polarization 
are still found. Under this interpretation, the possible RM 
gradient found would be parallel to the local direction of 
the jet. This could still be explained in terms of a shock, if 
we assume that we have a strong compression shock where 
the E-field is enhanced, followed by a dissipation zone pro- 
ducing the RM gradient. This shock might have been form 
recently, based on the proper motion of this component, 
w hich may expla i n wh y such a gradient was not observed 
m iHovatta et all (|2012l i as their observations were done in 
September 2006, more than two years prior to ours. 



However, we are currently una ble to distingu i sh be - 
tween these scenarios. According to iHovatta et al.l (|2009l l, 
the viewing angle for this source is only 2.5°, and so we are 
looking at it almost head-on, which makes it difficult to ad- 
equately study its de-projected features. Also, observations 
with more sensitivity would be necessary at higher frequen- 
cies to find out the actual direction of the jet at different 
scales. Further observations of the knot studied here will 
be necessary in order to understand its motion and E-field 
evolution in the forthcoming years. 



5 CONCLUSIONS 

We have obtained polarimetric VLBA observations of 8 
AGNs and presented their rotation measure analysis. Fara- 
day rotation is found in the VLBA cores with typical values 
for the rotation measure on the order of some thousands 
of rad/m^, except for the sources 1633+382 and 1954+513, 
with core RMs up to 22 krad/m^ and -13 krad/m^. Due to 
sensitivity limits, jet RM is detected only in the jets of two 
sources: 0906+430 and 1633+382. 

Core RMs found here are larger than in previous stud- 
ies performed in the same sources at lower frequencies. If 
we assume a d epend ence of the form RM oc zy", derived by 
Ijorstad et aD (|2007l ). we find an average of a = 3.6 ± 1.3, 
in agreement with estimations made for simil ar sources 
(lO'Sullivan k Gabuzdal[2009l : ljorstad et al.ll2007f ). However, 
most of the individual values are higher than the theoretical 
estimations that provide a = 2 under the assumption of a 
toroidal dominated magnetic field and electron density in 
equipartition. Different explanations for this include a time 
variability of the RM (supported in some sources by a sign 
change over different epochs) or a convolution effect causing 
RMs to be underestimated at lower frequencies. 

In order to discuss the reliability of the RMs found 
here we have performed simulations based on our multi- 
frequency data. Our simulations indicate that for RM vari- 
ations to be robust, they have to span over at least 1.5 beam 
sizes with a SNR> 3, although we have found that these con- 
ditions can be slightly loosened as stable gradients over a 
few beam sizes should not arise due to random variations or 
be am effects. These result s agr ee with previous disc ussions 
bv lTavlor Zavala! (|2010l ) and IHovatta et al.1 (l2012l ). 

The Faraday-corrected direction of the core magnetic 
field is, in general, either aligned or parallel with the direc- 
tion of the jet. This can be interpreted as an indication of an 
underlying helical magnetic field. Two sources seem not to 
follow this trend. In 1633+382 there were various possibili- 
ties including a change of the jet direction and in 1823+382 
we propose that the rotation of the polarization angle may 
be due to an unresolved component. This interpretation is 
supported by observations at higher frequencies and slices 
of the degree of polarization, revealing this component. 

In some sources (0133+476, 0745+241 and 0906+430) 
we observe two different regions for the polarization angle. 
Based on the previous interpretation, we suggest that this is 
due to the unaltered E-field in the core plus an additional 
component that, being still unresolved at 12 GHz, has moved 
away from the core enough to be detectable in polarization. 

We find that the RM gradient in the jet of 0906+430 
(and possibly 1633+382) is reliable. The combination of the 
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results from the i?M, the spine polarization structure and 
the concave shape of the degree of polarization across the 
jet seem to indicate the presence of a helical magnetic field 
in 0906+430 beyond 1 mas away from the core. The case 
for 1633+386 is unclear and different possibilities, such as a 
compression shock or diffusion of enhanced magnetic field, 
are given. 



APPENDIX A: RELIABILITY OF ROTATION 
MEASURE GRADIENTS 

The subject of rotation measure gradients in the last decade 
has been a controversial one. Due to the difficulty in address- 
ing the RM errors and the extend of the RM compared with 
the map beam sizes, the reliability of RM gradients is still 
debatable. Initially one might think that, when the size of 
the region where the observed rotation measure is compa- 
rable to the beam size, we cannot truly resolve the RM 
features and the gradients we find might be the spurious 
product of the integration of several random features within 
the restoring beam, and not real physical on es. In order to 
clarify the situation, iTavlor Zavalal (|2010l ) summarized a 
series of criteria in order to establish observational require- 
ments for a reliable RM gradient. We summarize and discuss 
these requirements here. 

As a gradient will exist whenever two values differ, it 
is clear that, at least, three independent measurements are 
needed. An initial (but too conservative) guess is that for 
each element of beam width one can only obtain a single 
independent value and hence, a gradient spanning along 
at least three beam widths is needed in order to speculate 
about a gradient. However, this requirement is too restric- 
tive. The position accuracy of polarized intensity images is 
below a small fraction of the synthesized beam size, because 
in the calibration process the same comp lex gains are ap - 
plied to both pa rallel and cross visibilities (|LepDanenlll995h . 
iLobanovl (l2005l ) argues that the minimum resolvable size is 
determined not only by the beam size but also by the signal- 
to-noise ratio, concluding that, in general, at SNR > 4, the 
minimum resolvable size is always smaller than the beam 
FWHM. For a typical SNR ^ 20, this implies a resolvable 
size of 1/3 of the beam size; i. e, ~ 0.25 mas for our case. 

iMurphv Gabuzdal (|2012l ) argue that an RM gradient 
can change its value when convolved with different beams, 
but the gradient itself is not destroyed even when the jet is 
unresolved, although they do no t take into account s puri- 
ous gradients arising from noise. iHovatta et aD (|2012l ) per- 
form a series of simulations specifically dealing with this 
issue that indicate that the gradient is already considered 
significant with a size of 1.5 beams, provided it exceeds 
3x RMS. This is consistent with similar criteria described 
bv lBroderick M cKinnev (2010) and in agreement with the 
discussion bv iLobanovl (|2005l ) . 

To study these aspects i n a more quantitative way we 
followed iHovatt a et al.l (| 20121 ) and performed a series of sim- 
ulations. In their work they used frequencies between 8 and 
15 GHz and explicitly pointed out that independent simula- 
tions should be done if a different set of frequencies were to 
be used, as it is the case here. Hence, we used their procedure 
for our case. 

First, the stokes I model of the CLEANed source was 



used. We then created fiducial Q and U models by setting 
them to be a known fraction of Stokes I (which leads to a 
constant EVPA and polarization across the source) . We then 
produced the UV data based on these models and added 
random independent noise with a normal distribution to the 
data. Using this method we produced 500 sets of UV data 
for each frequency, which were used then to derive new Q, 
U, polarized and respective noise maps. We checked that the 
resulting maps had noise levels and distributions statistically 
compatible with the real observations. 

We then created simulated RM and RM noise maps us- 
ing the values of the polarization angles from the simulated 
12,15 and 22 GHz maps, to which we added a random value 
between ±4° to account for EVPA calibration errors. Errors 
on the RM were obtained by the quadratic sum of differ- 
ent terms, namely the errors from the noise map, estimated 
from the variance- covariance matrix of the linear fit as the 
error of the slope, and the simulated propagated errors from 
the polarization maps. This is done because, although er- 
rors from the covariance matrix alone are adequate for each 
pixel, we still need to take into account that pixels are not 
independent. We note that we do not include the EVPA cal- 
ibration errors, as they will have no effect in RM gradients 
iMahmud etTahl (l2009bl \ which we are interested in. Figure 
Al shows a sample of a simulated RM map. 

Typical error in the RM is about 150 rad/m^, which 
is a bit smaller than the error in the real RM in our obser- 
vation at the same position 250 rad/m^). %Nonetheless, 
when we carefully check our estimations, we are not able to 
find a definitive cause for this. We estimate that there is a 
probability of ~5% to obtain a simulated RM with an error 
similar or larger than the observed one. Thus, it is possi- 
ble that our observation is within this case. It might also be 
possible that this difference arises due to underestimation of 
the non-thermal errors due to the procedure in creating the 
Q and U models, or the evaluation of the pixel dependence. 
After carefully checking out estimations we think this is not 
the case. 

Figure A2 shows the distribution of obtained RM values 
at a point in the jet for the source 1633+382 close to the 
area where the actual RM is observed. Values are about a 
third smaller than the ones in the actual observation, around 
few hundred rad/m^. The distribution of RM values found 
is as follows: the average and the most probable values are 
-8 rad/m^, and -9 rad/m^, respectively, consistent with a 
distribution that has no RM. The standard deviation of the 
distribution is 104 rad/m^. 

For each of the simulated RM maps, we took a slice 
across the jet, as shown in Figure Al by the straight line. 
Ideally, we would like to take different slices of various sizes 
separated by at least one beam width to ensure they are in- 
dependent. However, given the structure of the source, this is 
not possible, particularly for slices smaller than ^2.5 beams. 
In order to overcome this issue we have used the following 
procedure: 

We first set the longitude of a single slice to 1, 1.5, 2, 
2.5 and 3 beam sizes respectively (we note here that 1 beam 
corresponds to 0.75 mas), all centered around the jet spine. 
In total, we obtained 2500 slices that were used to calculate 
simulated RM gradients as follows: we fitted a simple line to 
the slice. We estimated the gradient error as the maximum 
between the error in the slope fit and the largest error bar in 
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Figure A5. Fraction of rotation measures along the slice in our 
simulation that exceeded 1 ("+" symbol, red in the online ver- 
sion), 2 ("x" symbol, green in the online version) and 3 (">|<" 
symbol, blue in the online version) standard deviations compared 
with size of the slice. 



Figure Al. Slices transverse to the jet in 1633+382 overlaid on 
a simulated RM map. Straight line indicates the slice used to 
obtain the simulated RM gradients. Dashed lines indicate slices 
that were taken to check for the consistency of the method. Color 
figure can be found in the online version of the journal. 
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Figure A2. Distribution of simulated rotation measures in the 
jet of 1633+382 



defined as above. Results can be seen in Figure A4, where 
the number of RM gradients are plotted in terms of their 
SNR. We note here that none of these RM gradients are 
real. However, when a slice with the size of only 1 beam 
is analyzed, there is still a significant number of spurious 
gradients with a =SNR> 3. The number of these gradients 
decreases dramatically as we study slices with 1.5 beam sizes 
or larger: only 4 sources for 1.5 beam sizes and none of them 
for larger ones had a > 3 in our simulation. 

This is clearly shown in Figure A5, where we plot the 
number of false positives (i.e, the proportion of sources that 
exceed la, 2a and 3a respectively) against the size of the 
slice. This fraction goes rapidly to zero as we increase the 
number of error limits a that are imposed. In general, la 
gradients will not be reliable even for large slices but a level 
of 2a will be enough when the width of the slices used is 
at least 2 beams. In most cases, even a 2a detection level 
over 1.5 beams will be reliable, with a chance of finding false 
gradients in only 0.8% of these. 

Summarizing, the results of our si mulations agree with 
the ones shown in 'Hova tta et aD (|2012l ) and the suggestions 
bv pTavlor & Zavala (201^ although we consider the latter 
too restrictive. While slices with 3 beam widths are optimal 
for a reliable detection of a RM gradient, in the most gen- 
eral case smaller slices can also provide a reliable gradient, 
provided the error treatment and SNR is adequate. 



the RM slice. Because this method does not provide inde- 
pendent measurements, we repeated this process for another 
two different slices shown by dashed lines in Figure Al as 
a consistency check. We find that the results are similar for 
all three slices. Figure A3 shows the distribution of RM gra- 
dients obtained by this method for the straight slice, which 
is taken along our observed RM. In all cases, the simulated 
RM is centered at ~ rad/m^, as expected, with the stan- 
dard deviation decreasing as we move to slices larger than 1 
beam size. 

To study when these RM gradients are significant, we 
have analyzed their SNR by dividing them by their error. 
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Figure A3. Distributions of the 2500 simulated rotation measure gradients. 
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Figure A4. Distributions of the 2500 simulated rotation measures gradients significance. 



REFERENCES 

Algaba, J. C, PhD Thesis, University College Cork (2010) 
Algaba, J. C, Gabuzda, D. C, Smith, P. S., MNRAS, 420, 
542 (2012) 

Algaba, J. C, et al. (2013), in preparation 

Aloy, M. A., Gomez, J. L., Ibanez, J. M., Marti, J. M., 



Miiller, E, ApJ 528, L85 (2000) 
Asada, K., Inoue, M., Uchida, Y., Kameno, S., Fujisawa, 

K., Iguchi, S., Mutoh, M., PASJ, 54, L39 (2002) 
Asada, K., Inoue, M., Kameno, S., Nagai, H., AJ, 675, 79 

(2008) 

Attridge, J. M.; Roberts, D. H.; Wardle, J. F. C, ApJ, 518, 
87 (1999) 



14 J. C. Algaba 



Attridge, J. M.; Wardle, J. F. C, Homan, D. C, ApJ, 633, 
85 (2005) 

Blandford, R. D. in Astrophysical Jets, ed. D. Burgarella, 
M. Livio, & C. P. O'Dea Cambridge: Cambridge Univ. 
Press (1993) 

Broderick, A. E., McKinney, J. C, AJ 725. 750 (2010) 
Contopoulos, I, Christodoulou, D. M., Kazanas, D. 

Gabuzda, D. C, ApJ 702, L148 (2009) 
Croke, S. M. k Gabuzda, D. C, 2008, MNRAS, 386, 619 
Gabuzda, D. C, Pushkarev, A. B, Cawthorne, T, V., MN- 
RAS 319, 1109 (2000) 
Gabuzda, D. C, Murray, E., Cronin, P., MNRAS 351, L89 
(2004) 

Gabuzda, D. C, in Proceedings for 'The 8th European 

VLBI Network Symposium' PoS(8thEVN)011 (2006) 
Gomez, J. L., Marscher, A. P., Alberdi, A, Jorstad, S., 

Agudo, I., NRAO VLBA Memo 30 (1992) 
Gomez, J. L., Roca-Sogorb, M., Agudo, I., Marscher, A. P., 

Jorstad, S. G., ApJ 733, 11 (2011) 
Hovatta, T., Valtaoja, T., Tornikoski, N,. Lahteenmaki, A., 

A&A, 494, 527 (2009) 
Hovatta, T., Lister, M. L., Aller, M. P., Aller, H. D., 

Homan, D. C, Kovalev, Y. Y., Pushkarev, A. B., 

Savolainen, T. , 2012 .arXiv: 1205.674 6 
Jorstad, S. G., Marscher, A. P., Stevens, J. A., Smith, P. S., 

Forster, J. R., Gear, W. K., Cawthorne, T. V., Lister, M. 

L., Stirhng, A. M., Gomez, J. L., Greaves, J. S., Robson, 

E. I., AJ, 134, 799 (2007) 

Kato, Y., Hayashi, M. R., Matsumoto, R., ApJ 600, 338 
(2004) 

Kudryavtseva, N. A., Britzen, S., Witzel, A.; Ros, E., 
Karouzos, M., Aller, M. F., Aller, H. D., Terasranta, H., 
Eckart, A., Zensus, J. A., A&A 526, 51 (2011) 

Laing, R. A., MNRAS 193, 439 (1980) 

Leppanen, K. J., Zensus, J. A., Diamond, P. J., AJ, 110, 
2479 (1995) 

Liu, Y.; Jiang, D. R.; Shen, Z.-Q.; Karouzos, M., A&A 522, 
A5 (2010) 

Lister, M. L., Homan, D. C, AJ 130, 1389 (2005) 

Lobanov, A. P., AfcA, 330, 79 (1998) 

Lobanov, A. P,, |astro-ph /0503225 (2005) 

Lovelace, R. V. E, Li, H., Koldoba, A. V., Ustyugova, G. 
v., Romanova, M. M., ApJ 572, 445 (2002) 

Lynden-Beh, D, Boily, C, MNRAS 267, 146 (1994) 

Lynden-Bell, D, MNRAS 279, 389L (1996) 

Mahmud, M, Gabuzda, D. C, in Proceedings for 'The 9th 
European VLBI Network Symposium on The role of VLBI 
in the Golden Age for Radio Astronomy and EVN Users 
Meeting' PoS(IX EVN Symposium)011 (2009a) 

Mahmud, M., Gabuzda, D. C, Bezrukovs, V., MNRAS, 
400, 349, 2-12 (2009b) 

Marscher, A. P., Jorstad, S. G., Larionov, V. M., Aller, 
M. F., Aller, H. D., Lahteenmaki, A., Agudo, I., Smith, 
P. S., Gurwell, M., Hagen- Thorn, V. A., Konstantinova, 
T. S., Larionova, E. G., Larionova, L. V., Melnichuk, 
D. A., Blinov, D. A., Kopatskaya, E. N., Troitsky, I. S., 
Tornikoski, M., Hovatta, T., Schmidt, G. D., D'Arcangelo, 

F. D., Bhattarai, D., Taylor, B., Olmstead, A. R., Manne- 
Nicholas, E., Roca-Sogorb, M., Gomez, J. L., McHardy, I. 
M., Kurtanidze, O., Nikolashvih, M. G., Kimeridze, G. N., 
Sigua, L. A.; ApJ 710, 126 (2010) 

McKinney, J. C, MNRAS 368, 1561 (2006) 



O'Sullivan, S. O. P., Gabuzda, D. C, MNRAS 393, 429 
(2009) 

Murphy, E., Gabuzda, D. C, JPhCS, 355, 012009 (2012) 
O'Sullivan, S. O. P., Gabuzda, D. C, ASPC, 402, 209 
(2009) 

Owen, F.N., Hardee, P. E., Cornwell, T. J., ApJ 340, 698 
(1989) 

Pacholczyk, A. G., Radio Astrophysics. Freeman & Co., 

San Francisco (1970) 
Papageorgiou, A. PhD Thesis, University of Central Lan- 

canshire (2006) 
Pushkarev, A. B., Gabuzda, D. C, Vetukhnovskaya Yu. N, 

Yakimov, V. E., MNRAS 356, 859 (2005) 
Roberts, D. H., Wardle, J. F. C, ApJ, 759, 35 (2012) 
Uchida Yu., Kigure H., Hirose S., Nakamura M. & Cameron 

R., ApJ, 600, 88 (2004) 
Taylor, G. B. & Zavala, R., ApJL, 722, 183 (2010) 
Venturi, T, Taylor, G. B., AJ 118, 1931 (1999) 
Zavala, R.T. & Taylor, G. B., ApJ, 550, 147 (2001) 
Zavala, R.T. & Taylor, G. B., ApJ, 589, 126 (2003) 
Zavala, R.T. & Taylor, G. B., ApJ, 612, 749 (2004) 
Zavala, R.T. & Taylor, G. B., ApJ, 626, 73 (2005) 



